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Abstract.  Transport  properties  of  impurities  immersed  in  a  granular  gas  undergoing  homogenous  cooling  state  are  studied. 
The  results  are  obtained  from  the  Boltzmann-Lorentz  kinetic  equation  for  inelastic  Maxwell  models  in  d  dimensions.  The 
kinetic  equation  is  solved  by  means  of  the  Chapman-Enskog  method  up  to  first  order  in  the  density  gradient  of  impurities. 
The  mass  and  heat  fluxes  associated  with  impurities  are  determined  and  the  corresponding  transport  coefficients  are  identified. 
Both  transport  coefficients  (diffusion  and  Dufour  coefficients)  are  exactly  obtained  in  terms  of  the  coefficients  of  restitution 
for  the  impurity-gas  and  gas-gas  collisions  as  well  as  the  ratios  of  mass  and  diameters.  The  results  are  compared  with  those 
obtained  for  inelastic  hard  spheres  in  the  leading  Sonine  approximation  and  by  means  of  Monte  Carlo  simulations.  The 
comparison  shows  good  agreement  between  both  interaction  models,  especially  in  the  case  of  the  diffusion  coefficient. 


1.  INTRODUCTION 

Granular  fluids  are  usually  modelled  by  an  idealized  system  of  smooth  hard  spheres  with  inelastic  collisions.  In  the 
low-density  limit,  the  Boltzmann  equation  has  been  conveniently  modified  to  account  for  the  inelasticity  of  binary 
collisions  and  the  Navier-Stokes  transport  coefficients  have  been  computed  by  solving  the  corresponding  kinetic 
equation  by  means  of  the  Chapman-Enskog  method  [1].  This  task  has  been  carried  out  in  the  past  few  years  in  the 
case  of  a  monocomponent  gas  [2]  as  well  as  for  multicomponent  systems  [3].  In  both  cases,  as  happens  for  elastic 
collisions,  all  the  transport  coefficients  are  given  in  terms  of  the  solutions  to  linear  integral  equations  which  are  solved 
approximately  by  taking  the  leading  terms  in  a  Sonine  polynomial  expansion.  In  spite  of  this  approximation,  the  kinetic 
theory  results  compare  quite  well  with  those  obtained  from  numerical  solutions  of  the  Boltzmann  equation  by  using 
the  direct  simulation  Monte  Carlo  (DSMC)  method  [4].  The  fact  that  the  agreement  between  kinetic  theory  and  DSMC 
results  even  occurs  for  strong  dissipation  is  a  further  testimony  of  the  utility  of  the  kinetic  theory  and  the  hydrodynamic 
descriptions  for  granular  fluids  beyond  the  weak  dissipation  limit. 

The  main  mathematical  difficulty  in  solving  the  Boltzmann  equation  for  inelastic  hard  spheres  (IHS)  comes  from 
the  form  of  the  collision  rate,  which  is  proportional  to  the  magnitude  of  the  relative  velocity  of  the  two  colliding 
particles.  In  the  case  of  elastic  collisions,  an  alternative  to  overcome  this  problem  is  to  consider  the  repulsive  Maxwell 
potential,  for  which  the  collision  rate  is  independent  of  the  relative  velocity.  This  property  simplifies  the  mathematical 
structure  of  the  Boltzmann  collision  operator  and  allows  one  to  get  exact  solutions  of  the  Boltzmann  equation  for  far 
from  equilibrium  situations  [5].  In  the  case  of  inelastic  gases,  the  so-called  inelastic  Maxwell  models  (IMM)  [6]  have 
been  widely  used  in  the  past  few  years  as  a  toy  model  to  characterize  the  effect  of  the  inelasticity  of  collisions  on 
the  properties  of  granular  gases.  The  IMM  share  with  elastic  Maxwell  molecules  the  property  that  the  collision  rate  is 
velocity  independent  but  their  scattering  rules  are  the  same  as  for  IHS. 

Most  of  the  studies  performed  in  the  context  of  IMM  have  been  devoted  to  homogeneous  states,  mainly  to  the 
analysis  of  the  overpopulated  high  energy  tails  [7].  Much  less  is  known  for  inhomogeneous  situations,  especially  for 
multicomponent  systems.  One  of  the  few  studies  considering  mixtures  of  IMM  has  been  carried  out  by  the  author  [8] 
in  the  uniform  shear  flow  problem.  The  rheological  properties  (shear  and  normal  stresses)  of  the  mixture  have  been 
exactly  evaluated  in  terms  of  dissipation  and  parameters  of  the  mixture  and  surprisingly,  the  results  of  IMM  show  an 
excellent  agreement  with  those  obtained  for  IHS  by  means  of  Grad’s  approximation  and  Monte  Carlo  simulations. 
This  fact  stimulates  the  search  for  exact  solutions  for  IMM  which  can  be  confronted  with  the  results  obtained  for  IHS 
by  using  approximate  analytical  methods  and  computer  simulations. 

The  goal  of  this  paper  is  to  determine  the  diffusion  coefficient  of  impurities  immersed  in  a  gas  of  inelastic  Maxwell 
molecules.  The  gas  is  in  the  homogeneous  cooling  state  (HCS),  namely,  a  homogeneous  state  where  the  temperature 
uniformly  decreases  in  time  due  to  collisional  cooling.  Diffusion  is  generated  in  the  system  by  the  presence  of  a 
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weak  concentration  gradient,  which  also  induces  a  heat  flux.  Both  mass  and  heat  fluxes  define  the  relevant  transport 
coefficients  of  the  problem:  the  diffusion  coefficient  and  the  Dufour  coefficient.  Here,  I  get  explicit  expressions  for 
these  coefficients  by  solving  the  corresponding  Boltzmann-Lorentz  equation  for  impurities  by  means  of  the  Chapman- 
Enskog  method  through  the  first  order  in  the  concentration  gradient.  These  transport  coefficients  are  given  in  terms 
of  the  coefficients  of  restitution  for  the  impurity-gas  and  gas-gas  collisions  as  well  as  the  masses  and  diameters  of 
impurities  and  gas  particles.  The  dependence  of  these  coefficients  on  the  parameters  of  the  system  is  illustrated  and 
compared  with  known  results  derived  for  IHS  [3,  9], 


2.  DESCRIPTION  OF  THE  PROBLEM 

Let  us  consider  a  dilute  gas  of  inelastic  Maxwell  molecules  of  mass  m,  diameter  cr,  and  interparticle  coefficient 
of  restitution  a.  In  the  low-density  regime,  its  velocity  distribution  function  /(v)  obeys  the  nonlinear  Boltzmann 
equation  (BE).  We  assume  that  the  gas  is  in  the  HCS,  namely,  a  homogeneous  solution  of  the  BE  in  which  all  the  time 
dependence  of  /(v,t)  is  through  the  granular  temperature  T(t).  In  this  case,  the  BE  can  be  written  as 

(l) 

where  the  Boltzmann  collision  operator  J[y\f,f]  for  IMM  is 

J[vi\f,f]  =  J  dy2  J  do  [a_1/(vi,f)/(v2,f)-/(vi,f)/(v2,0]  •  (2) 

Here,  n  is  the  number  density  of  the  gas  particles,  co  is  an  effective  collision  frequency  (to  be  chosen  later),  £lrj  = 
2n,U2/V(d/2)  is  the  total  solid  angle  in  d  dimensions,  and  a  is  a  unit  vector  directed  along  the  centers  of  the 
two  colliding  spheres.  Upon  writing  Eq.  ( 1 ),  I  have  taken  into  account  the  balance  equation  for  the  temperature 
d,T  =  —  £,  where  the  cooling  rate  C  is  given  by  [10] 

C(a)  =  ^fm.  (3) 

The  explicit  form  of  the  velocity  distribution  function  /(v,  t )  is  not  explicitly  known,  although  its  velocity  moments 
can  be  explicitly  computed  [10]. 

We  assume  now  that  some  impurities  of  mass  «?o  and  diameter  Oq  are  added  to  the  system.  Given  that  its  molar 
fraction  is  negligible,  the  state  of  the  granular  gas  is  not  disturbed  by  the  presence  of  impurities  and  so  the  velocity 
distribution  function  /(v,  t )  obeys  the  BE  (1).  Furthermore,  collisions  among  impurities  themselves  can  neglected  ver¬ 
sus  the  impurity-gas  collisions,  which  are  characterized  by  a  coefficient  of  restitution  do.  As  said  in  the  Introduction, 
I  want  to  determine  the  transport  properties  of  impurities  generated  in  the  system  by  the  presence  of  a  weak  concen¬ 
tration  gradient  Vxo,  xo  =  no/n  being  the  mole  fraction  of  impurities.  Under  these  conditions,  the  velocity  distribution 
function  fo(r. \.t)  of  impurities  verifies  the  Boltzmann-Lorentz  equation 

dtfo+y-  V/0  =J[v\fo,f],  (4) 

where 

•/[vi|/o,/]  =  J dy2  J d o  [a()~ 1  /o (r , V; ; f )/(r , V, ; t )  —  /o (r , vi ; f )/(r, V2 ; f )]  ,  (5) 

coo  is  an  effective  collision  frequency  (to  be  chosen  also  later)  for  impurity-gas  collisions,  and 

v'i  =  v i - - — (1  +  «Q  ')  (<J-gi2)<7,  v'2=v2H - (l  +  a^1)  (CT-gi2)CT  ,  (6) 

m  +  mo  m  +  mo 

where  gi2  =  vj  —  v2  is  the  relative  velocity  of  the  colliding  pair.  The  partial  temperature  of  impurities  Tq  is  defined  as 

2 no To  =  J  dv”^v2/o(v),  (7) 

where  no  is  the  number  density  of  impurities.  In  contrast  to  what  happens  for  normal  fluids,  if  impurities  and  gas 
particles  are  mechanically  different,  the  partial  temperature  To  can  differ  from  the  granular  temperature  T  when 
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dissipation  is  present.  The  failure  of  energy  equipartition  in  granular  fluids  has  even  observed  in  real  experiments 
[1 1].  Momentum  and  energy  are  not  collisional  invariants  of  the  Boltzmann-Lorentz  collision  operator  J[fo,f\ ■  Only 
the  number  density  of  impurities  is  conserved: 


dtn0  + 


V-jo 

mo 


=  0, 


(8) 


where  the  mass  flux  jo  is  defined  as 

jo  =  J  d\m0yfo{y)-  (9) 

In  the  absence  of  diffusion,  impurities  are  also  in  HCS  and  Eq.  (4)  becomes 

-(v/o)  =J[y\fo,f],  (10) 

where  use  has  been  made  of  the  relation  T()  1  dt  To  =  —  £0,  where  Co  is  the  cooling  rate  associated  with  the  partial 
temperature  To  of  impurities.  Its  expression  is  [8] 

Co  =  — M(1  +<*0)  [l  —  ^(1  +  ofo)(l  +  0)j  >  (11) 


where  fi  =  m/(m  +  mo )  and  Q  =  moT /mTo  is  the  mean  square  velocity  of  gas  particles  relative  to  that  of  impurities. 
The  balance  equations  for  the  temperatures  T  and  To  lead  to  the  following  time  evolution  equation  for  the  temperature 
ratio  y  =  Tq/T: 

Yl^tY=  C  —  Co-  (12) 


The  fact  that  fo  depends  on  time  only  through  T  ( t )  necessarily  implies  that  the  temperature  ratio  y  must  be  independent 
of  time,  and  so  Eq.  (12)  gives  the  HCS  condition  C (f)  =  Co(0-  The  impurity  equilibrates  to  a  common  HCS  with 
different  temperatures  for  the  impurity  and  gas  particles.  This  implies  a  breakdown  of  the  energy  equipartition. 

In  order  to  get  the  explicit  dependence  of  y  on  the  parameter  space,  one  still  needs  to  fix  the  parameters  co  and 
C0q.  To  optimize  the  agreement  with  the  IHS  results,  the  cooling  rates  C  and  Co  of  IMM  given  by  Eqs.  (3)  and  (11), 
respectively,  are  adjusted  to  be  the  same  as  the  ones  given  for  IHS  [12],  Since  C  and  Co  are  not  exactly  known  for  IHS, 
one  can  estimate  them  by  taking  their  local  equilibrium  approximations.  In  this  case,  the  collision  frequencies  co  and 
C0o  are  given  by 


o.d  (  1+0 

Vet/  \ 


(13) 


where  a  =  (cr  +  cto) / 2  and  v  =  na'1  1  \/2T / m  is  an  effective  collision  frequency.  With  these  choices,  the  dependence 
of  the  temperature  ratio  Tq/T  on  the  parameters  of  the  system  given  for  IMM  presents  an  excellent  agreement  with  the 
one  found  for  IHS,  even  for  strong  dissipation  [13]. 


3.  TRANSPORT  PROPERTIES  OF  IMPURITIES 


The  aim  of  this  Section  is  to  get  the  transport  properties  (mass  and  heat  fluxes)  of  impurities  in  the  limit  of  small 
concentration  gradient  Vjt'o-  These  fluxes  are  determined  by  solving  the  Boltzmann-Lorentz  equation  (4)  from  the 
Chapman-Enskog  method  [1]  adapted  to  dissipative  systems.  Thus,  I  look  for  a  normal  solution  in  which  all  the 
space  and  time  dependence  of  fo  is  through  the  hydrodynamic  fields.  In  our  problem,  this  means  that  fo(r. \.r)  = 
/o[v|jt'o(f),  T(t)\.  The  normal  solution  is  generated  by  expanding  fo  in  powers  of  V.v(l  as 

/o  =  /0(0)+£/<51)  +  ---,  (I4) 

where  each  factor  e  corresponds  to  the  implicit  factor  Vxo-  The  time  derivative  is  also  expanded  as  d,  =  r),'  0  :  +  e  d,  1 + 
■  ■  ■ ,  where  the  action  of  the  different  operators  c)\ k :  can  be  obtained  from  the  balance  equation  (8).  They  are  given  by 


4<0)-*o  =  0, 


Y7  •(*-!) 

d,(0)T  =  -T £,  dt(k)xo  =  --  ’Jo 


moil 


dt{k)T  =  0,  k>  1, 


(15) 
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FIGURE  1.  Plot  of  the  reduced  diffusion  coefficient  D(a)/D(  1)  as  a  function  of  the  coefficient  of  restitution  a  =  a o  in  the 
three-dimensional  case  for  Gq/ a  =  2  and  mo/m  =  8.  The  solid  line  corresponds  to  the  exact  result  obtained  here  for  IMM  while 
the  dashed  line  is  the  result  derived  for  IHS  in  the  second  Sonine  approximation.The  symbols  refer  to  Monte  Carlo  simulations  for 
IHS. 


where 


jo°  =  J 


(16) 


The  zeroth-order  approximation  /q  is  the  solution  of  Eq.  (10)  but  taking  into  account  now  the  local  dependence  on 
the  mole  fraction  xo .  Since  /q  is  isotropic,  it  follows  that  the  flux  of  impurities  vanishes  at  this  order,  i.e.,  j0°;  =  0, 
and  so  d} 1  xo  =  0.  To  first  order  in  e,  one  has  the  kinetic  equation 


(0) 


-(J-fV) 


\dxo“ 


fo  vi  '  Vx-o. 


(17) 


The  second  equality  follows  from  the  balance  equations  (15)  and  the  space  dependence  of  /q  through  xo-  The  mass 

flux  jg  and  the  contribution  to  the  heat  flux  coming  from  impurities  qj,1 :  can  be  easily  obtained  from  Eq.  (17).  The 
heat  flux  is  defined  as 

qo1)  =  /rfv'?v2v/o(1)-  (18) 

To  get  these  fluxes,  one  multiplies  both  sides  of  Eq.  (17)  by  niov  and  !, niov2v  and  integrates  over  v.  The  result  for 
the  mass  and  heat  fluxes  is,  respectively 


(a,™  +  v„)q™  _ 

where  p  =  nT  is  the  hydrostatic  pressure  and 


(d/0)  +  vd)  j^1}  =  - ypVx0 , 

(l  +  f)y2Vx0, 


d  pT 
d  +  2  mo 


c  o  = 


d(d  +  2)  [4 7i()T02 


J  d\v4fo(\)  — 


d(d  2) 


Upon  writing  Eqs.  (19)  and  (20),  use  has  been  made  of  the  collisional  moments  [13] 

j  d\mo\J[f^\f}  =  -vj^, 


dy  ^v2v7[/q1},/]  =  -v9q(l/;  -  pi{01 


,(i) 


(19) 

(20) 

(21) 

(22) 

(23) 
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FIGURE  2.  Plot  of  the  reduced  Dufour  coefficient  Lq(a)  as  a  function  of  the  coefficient  of  restitution  a  =  <Zo  in  the  three- 
dimensional  case  for  <7o  =  <7,  and  two  different  values  of  the  mass  ratio:  mo/m  =  2  and  mo/m  =  8.  The  solid  lines  correspond  to 
the  exact  results  obtained  here  for  IMM  while  the  dashed  lines  are  the  results  derived  for  IHS  in  the  first  Sonine  approximation. 


where 


Vd  = 


s/nd  \a 


o'd~l 


1  +  0 
0 


1/2 


ya  =  - 


Vd 

d  - b  2 


+  Oo)  [d  +  8  —  3/i(  1  H-  cxq )]  —  ?>{d  +  2) } , 

VdT 


P  —  -z — iz(l  «b)  \d  —  3/i(l  +  ao)  +  2] . 

Zm 


The  mass  and  heat  fluxes  have  the  forms 


j(l')  =  -moDVx  o, 

,(U 


(24) 

(25) 

(26) 

(27) 

(28) 


q0' 1  =  -T2DgVx o, 

where  D  is  the  diffusion  coefficient  and  Dq  is  a  cross  transport  coefficient  measuring  energy  transport  due  to  diffusion 
of  impurities.  Dimensional  analysis  requires  that  D  °c  7 1  /2  and  Dq  °c  7'  1  2.  Consequently, 


=  -QTdT jW  =  ^DVx o,  -  +:C7»/  qi,l!  =  \^T2DqWx0. 


(29) 


Inserting  Eqs.  (29)  into  Eqs.  (19)  and  (20),  respectively,  the  transport  coefficients  D  and  Dq  can  be  easily  identified 
from  (27)  and  (28).  The  result  is 

D=— - 7 ,  ..  ,  (30) 


Dn  = 


m0v  y*  -  ’ 

n  (mo/m)p*D* -^(l  +  cf)y2 
niQ  V  f£*-V* 


(31) 


2^  yq 

Here,  =  vrf/v,  /3*  =  /3 m/(Tv),  v*  =  v9/v,  and  £*  =  £/v. 

A  useful  way  of  characterizing  the  heat  transport  in  a  mixture  due  to  a  concentration  gradient  is  through  the  Dufour 
coefficient.  In  our  problem,  the  (dimensionless)  Dufour  coefficient  Lq  is  defined  by  the  relation 


T  ft)  d  +  2i^  »r2.v. 

.  J , ,  —  q()  -  T  —  I.q  \  \(  i . 

*  U  2  Mi()  7H()V  1 

The  coefficient  Lq  can  be  easily  written  in  terms  of  the  coefficients  D  and  Dq  as 


Lq  — 


wiqV 


T°q-d-/D 


(32) 


(33) 
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In  the  elastic  case  (a  =  Oq  =  1),  £  =  0,  Tq  =  T,  and  as  expected  Lq  =  0.  Figures  1  and  2  show  the  dependence  of 
the  reduced  coefficients  D(n)  j D(  I )  and  Lq(a),  respectively,  on  the  (common)  coefficient  of  restitution  a  =  Oq  for 
different  systems.  Here,  D(  1 )  is  the  elastic  value  of  the  diffusion  coefficient.  In  the  case  of  the  diffusion  coefficient,  it  is 
apparent  that  the  results  for  IMM  and  IHS  are  practically  indistinguishable  over  the  whole  range  of  values  of  a  studied. 
Moreover,  both  theories  present  an  excellent  agreement  with  Monte  Carlo  simulations.  More  discrepancies  between 
IMM  and  IHS  are  found  in  the  case  of  the  Dufour  coefficient  Lq,  especially  with  increasing  dissipation.  However,  at  a 
qualitative  level,  the  IMM  predictions  compare  well  with  the  ones  derived  for  IHS. 


4.  CONCLUDING  REMARKS 

Diffusion  of  impurities  in  a  dilute  granular  gas  of  IMM  undergoing  HCS  has  been  analyzed.  The  relevant  transport 
coefficients  of  the  problem  (diffusion  and  Dufour  coefficients)  have  been  determined  from  a  Chapman-Enskog  solution 
[1]  of  the  Boltzmann-Lorentz  equation  up  to  first  order  in  the  concentration  gradient.  The  expressions  (30),  (31),  and 
(33)  derived  here  for  the  diffusion  D  and  Dufour  Lq  coefficients  are  exact  (within  the  context  of  IMM)  and  constitute 
the  main  goal  of  this  paper.  This  contrasts  with  the  previous  results  derived  for  IHS  [3,  9],  where  the  coefficients  D 
and  Lq  were  obtained  by  considering  the  leading  terms  in  a  Sonine  polynomial  expansion  of  the  velocity  distribution 
function  of  impurities.  To  make  contact  with  the  results  reported  for  IHS,  the  collision  frequencies  co  and  C0q  (which 
can  be  seen  as  free  parameters  of  the  model)  need  to  be  chosen.  Here,  I  have  taken  these  frequencies  to  reproduce  the 
cooling  rates  C  and  £o  of  IHS  in  the  local  equilibrium  approximation.  With  this  choice,  the  comparison  with  known 
results  for  IHS  shows  that  in  general,  the  IMM  predictions  are  reasonably  good  for  not  too  large  dissipation,  especially 
for  the  diffusion  coefficient  which  is  related  to  a  first-degree  velocity  moment  of  the  velocity  distribution  function. 
In  the  case  of  the  Dufour  coefficient,  although  the  IMM  capture  qualitatively  well  the  trends  obtained  for  IHS,  the 
disagreement  between  both  interaction  potentials  is  more  significant  than  the  one  found  for  the  diffusion  coefficient 
D.  In  this  sense,  more  discrepancies  between  IMM  and  IHS  are  expected  when  one  considers  higher  degree  velocity 
moments. 
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